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This paper considers magnitude, asymptotics and duration of drawdowns for some Levy pro¬ 
cesses. First, we revisit some existing results on the magnitude of drawdowns for spectrally 
negative Levy processes using an approximation approach. For any spectrally negative Levy 
process whose scale functions are well-behaved at O-f, we then study the asymptotics of draw¬ 
down quantities when the threshold of drawdown magnitude approaches zero. We also show that 
such asymptotics is robust to perturbations of additional positive compound Poisson jumps. Fi¬ 
nally, thanks to the asymptotic results and some recent works on the running maximum of Levy 
processes, we derive the law of duration of drawdowns for a large class of Levy processes (with 
a general spectrally negative part plus a positive compound Poisson structure). The duration 
of drawdowns is also known as the “Time to Recover” (TTR) the historical maximum, which 
is a widely used performance measure in the fund management industry. We find that the law 
of duration of drawdowns qualitatively depends on the path type of the spectrally negative 
component of the underlying Levy process. 
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1. Introduction 

Drawdowns relate to an investor’s sustained loss from a market peak. It is one of the most 
frequently quoted indices for downside risks in the fund management industry. Drawdown 
quantities appear in performance measures such as the Calmar ratio, the Sterling ratio, 
the Burke ratio, and others; see, for example, Schuhmacher and Eling [34] for a collection 
of such drawdown-based performance measures. Eurthermore, drawdown problems have 
drawn considerable theoretical and practical interest in various research areas including 
probability, finance, risk management, and statistics; see Section 1.1 for a brief literature 
review. 
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In this paper, we consider a one-dimensional Levy process X = {Xt,t > 0} defined on 
, t ^ a filtered probability space satisfying the nsual conditions. The 

drawdown process of X is defined as 

Yt = Mt-Xt, t>0, 

where Mt = supg<„<( X^ is the running maximum (historical peak) of X at time t. Let 

Ta = inf{t > 0 : Ft > a}, 

be the first time the magnitude of drawdowns exceeds a pre-specified threshold a > 0. 
Given that (maxo<s<t Yg > a) = (tq < t) P-a.s., the distributional study of the maximum 
drawdown in magnitude is equivalent to the study of the stopping time Tq. 

However, from a risk management standpoint, the magnitude itself is not sufficient 
to provide a comprehensive risk evaluation of extreme drawdown risks. For instance, for 
extreme risks such as tornado and flooding, it is natural to also investigate the frequency 
and the duration of drawdowns. Landriault et al. [22] recently studied the frequency 
of drawdowns for a Brownian motion process by defining two types of drawdown time 
sequences depending on whether a historical running maximum is reset or not. In this 
paper, we will consider the duration of drawdowns, also known as “Time to Recover” 
(TTR) the historic running maximum in the fund management industry. 

Mathematically, the duration of drawdowns of a stochastic process can be considered 
as the length of excursions from its running maximum. For t > 0, let Gt '■ = sup{0 < s <t: 
Yg = 0} be the last time the process Y is at level 0 (or equivalently X = M) before or at 
time t. The drawdown duration at time t is therefore t — Gt- We then define a stopping 
time 


? 7 f, = inf{t > 6 : t — Gt > 6 }, ( 1 . 1 ) 

that is the first time the duration of drawdowns exceeds a pre-specified time threshold 
b > 0. Equivalently, the event (fjt, > t) implies that the maximum duration of drawdowns 
before time t is shorter than b. 

The stopping time rjb is related to the so-called Parisian time, which is the first time 
the length of excursions from a fixed spatial level (rather than its running maximum) ex¬ 
ceeds a pre-specified time threshold; see, for example, Chesney et al. [9] and Czarna and 
Palmowski [11]. Further, Loeffen et al. [26] provided an unified proof to derive the prob¬ 
ability that the Parisian time occurs in an infinite time horizon (known as the Parisian 
ruin probability in actuarial science) for spectrally negative Levy processes. Notice that, 
in contrast to the Parisian time, the stopping time r]i, is almost surely finite (e.g., page 
105 of Bertoin [4]), which motivates us to study the Laplace transform (LT) of rg, in 
this paper. Another related concept is the so-called red period of the insurance surplus 
process; see Kyprianou and Palmowski [21]. The red period corresponds to the length 
of time an insurance surplus process shall take to recover its deficit at ruin. But it is 
different than the distributional study of 77 b, especially when X has no negative jumps 
(e.g.. Brownian motion). 
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1.1. Literature review on drawdowns 

Taylor [36] first derived the joint Laplace transform of and for Brownian motion 
processes. Later on, it was generalized by Lehoczky [24] to time-homogeneous diffusion 
processes. Douady et al. [13] and Magdon et al. [28] derived infinite series expansions 
for the distribution of Ta for a standard and drifted Brownian motion, respectively. For 
spectrally negative Levy processes, Mijatovic and Pistorius [29] obtained a general sex¬ 
tuple formula for the joint Laplace transform of Ta and the last passage time at level Mr,, 
prior to Ta, together with the joint distribution of the running maximum, the running 
minimum, and the overshoot of Y at Ta- Also, an extensive body of literature exists on 
the dual of drawdowns, drawups, which measure the increase in value of an underlying 
process from its running minimum; see, for instance, Pistorius [30], Hadjiliadis and Vecef 
[16], Pospisil et al. [32], and Zhang and Hadjiliadis [38, 39]. 

In finance and risk management, researchers have devoted considerable effort in as¬ 
sessing, managing, and reducing drawdown risks. For instance, Grossman and Zhou [15] 
examined a portfolio selection problem subject to drawdown constraints. Cvitanic and 
Karatzas [10] extended the discussion to multiple assets. Chekhlov et al. [8] proposed a 
new family of risk measures and studied problems of parameter selection and portfolio 
optimization under the new measures. Pospisil and Vecer [31] invented a new class of 
Greeks to examine the sensitivity of investment portfolios to drawdowns. Garr et al. [5] 
designed some European-style digital drawdown insurance contracts and proposed semi¬ 
static hedging strategies using barrier options and vanilla options. Other recent works 
on drawdown insurance are Zhang et al. [37, 40], among others. 

In addition, a few a priori unrelated problems in finance and insurance are also closely 
connected to the drawdown problematic. For instance, the pricing of Russian options 
(e.g., Shepp and Shiryaev [35], Asmussen et al. [1] and Avram et al. [2]), and the op¬ 
timal dividend models with “reflecting barriers” (e.g., Avram et al. [3], Kyprianou and 
Palmowski [21], and Loeffen [27]) are two common examples. 

1.2. Objective and structure 

In this article, we begin by developing an approximation technique in the spirit of 
Lehoczky [24] to revisit several known LT results on the magnitude of drawdowns of 
spectrally negative Levy processes via basic fluctuation identities. 

Second, as the threshold of drawdown magnitude a j, 0, we examine the asymptotic 
behavior of those LTs for any spectrally negative Levy process whose scale functions are 
well-behaved at 0-1- (see Assumption 4.1 below). We also show that such asymptotics are 
robust with respect to the perturbation of arbitrary positive compound Poisson jumps, 
and hence obtain the asymptotics of drawdown estimates for a class of Levy models with 
two-sided jumps. 

Finally, we study the duration of drawdowns via the LT of rj},. First, an approximate 
scheme for the LT of ry, is developed. To obtain a well-defined limit, we turn our problem 
to the behavior of the density of the running maximum process M and the convergence 
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of some potential measure of the drawdown process Y. Thanks to the asymptotic results 
obtained and some recent works on the distribution of the running maxima of Levy 
processes (e.g., Chaumont [6], Chaumont and Malecki [7], and Kwasnicki et al. [19]), we 
obtain the law of rih in terms of the right tail of the ascending ladder time process for 
a class of Levy process with two-sided jumps (a general spectrally negative part plus a 
positive compound Poisson structure). 

The rest of the paper is organized as follows. In Section 2, we review the scale function 
of spectrally negative Levy processes and the ascending ladder process of a general Levy 
process. In Section 3, we revisit some known LT results on the magnitude of drawdowns 
of spectrally negative Levy processes based on an approximation approach. The asymp¬ 
totic behavior of these LTs for small threshold is studied in Section 4, where we also 
examine the asymptotic behavior in the presence of positive compound Poisson jumps. 
In Section 5, the LT of rji, is derived for a large class of Levy processes with two-sided 
jumps. Some explicit examples are presented in Section 6. For completeness, some results 
on the extended continuity theorem are presented in the Appendix. 


2. Preliminaries 

In this section, we briefly introduce some preliminary results for Levy processes. Readers 
are referred to Bertoin [4] and Kyprianou [20] for a more detailed background. 

For ease of notation, throughout the paper, we let R = (— 00 , 00 ), 1R+ = [0,oo) and 
1HI+ = {s G C : Re(s) > 0}. We denote by Pa, the law of a Levy process with Xq = a; G M. 
For brevity, we write P = Pq. The minimum of real numbers u, v is denoted by u Av = 
min{u,t;}. For a function /(•) on (0,oo) and xq G [0,oo], we write f{x) = o{g{x)) as 
X ^ xq for a positive function g{-) if lima;_>. 3 ;|, f{x)/g{x) = 0. 

2.1. Spectrally negative Levy processes and scale functions 

Consider a spectrally negative Levy process X = {Xt,t > 0}. Throughout the paper, we 
assume that jAlj is not a subordinator and hence 0 is regular for (0,oo) (see Definition 6.4 
and Theorem 6.5 of Kyprianou [20] for the definition and equivalent characterizations of 
the regularity). The Laplace exponent of X is given by 

•0(s) :=i log]E[e®^*] =-/rs-k -k / (e""^ - 1 - sxl{ 3 ,>_i})n(da;), (2.1) 

^ ^ J ( —CXD,0) 

for every s G H+. Here, ct > 0 and the Levy measure n(da;) is supported on (— 00 ,0) with 
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It is known that X has paths of bounded variation if and only if |a;|n(da;) < oo 

and cr = 0. In this case, we can rewrite (2.1) as 


1 ^( 3 ) = sd+ f (e'*® — I)n(dx), s > 0, 
( — 00 , 0 ) 


( 2 . 2 ) 


where the drift d := —/r + |a;|n(da;) > 0 as \X\ is not a subordinator. For any given 

g > 0, the equation ip{s) = q has at least one positive solution, and we denote the largest 
one by $(g). 

It is well known that t > 0} is a martingale for any c > 0. This gives rise 

to the change of measure 


dP'^ 

dP 




— QcXt-ip[c)t 


t>Q. 


(2.3) 


Under the new measure P^, X is still a spectrally negative Levy process, and its Laplace 
exponent is given by '0c(s) = '4’[s + c) — 'ip{c) for all s £ C such that s + c G 1HI+. 

For any g > 0, the g-scale function : K !->■ [0, 00 ) is the unique function supported 
on (0,oo) with Laplace transform 


/ e Ax = -, s > <i>(g). 

Ao.oo) V{s)-q 

It is known that is continuous and increasing on (0,oo). Henceforth, we assume 

that the jump measure n(da;) has no atom, then it follows that € C^(0,oo) (e.g.. 
Lemma 2.4 of Kuznetsov et al. [18]). Moreover, if the Gaussian coefficient cr > 0 then 
£ (17^(0,00) for all g > 0 (e.g.. Theorem 3.10 of Kuznetsov et al. [18]). The g-scale 
function is closely related to exit problems of the spectrally negative Levy process 
X with respect to first passage times of the form 

r+(-)=inf{t>0:Xt>(<)a:}, x £ M. 


Two well-known fluctuation identities of spectrally negative Levy processes are given 
below (e.g., Kyprianou [20], Theorem 8.1). For g > 0 and 0 < a; < a, we have 


E,[e-«^‘^ 1 


{T+<T-} 


W^'i\x) 
0 " “ IF(9)(a) 


(2.4) 


and 

]Ex[e 0 l{j,-<j,+}] = (2.5) 

where {x) = l+q (y) Ay. 

The following lemma gives the behavior of scale functions at 0+ and 00 ; see, for 
example, Lemmas 3.1, 3.2 of Kuznetsov et al. [18]. Relation (2.6) is from (3.13) of Egami 
et al. [14]. 
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Lemma 2.1. For any q>0, 


VF('?)(0+) 


0, if(T> 


1 

d' 

2 


if a > 0 or |a;|n(da;) = oo (unbounded variation), 

otherwise (bounded variation), 




oo, 

q + n(— oo, 0) 


if(T>0, 

if a = 0 and n(— oo, 0) = oo, 
if a = 0 and n(— oo, 0) < oo. 


and 



( 2 . 6 ) 


2.2. The ascending ladder process of general Levy processes 

In this subsection, we consider a general Levy process X = {Xt,t > 0} characterized by 
its characteristic exponent 

^'(s) := - jlogE[e*'*^*] =i/rs++ [ (1 - + isxl{| 3 ;|<i})n(da;), (2.7) 

c ^ Jb.\{o} 

for all s G K. If X has bounded variation, we can rewrite (2.7) as 

■^{s) = -isd+ [ (l-e*'’^)n(da;), (2.8) 

where the drift d := —/i — /o<| 3 ,|<i a;n(dx). 

The local time of X at its running maximum, denoted by L = {Lt,t > 0}, is a contin¬ 
uous, non-decreasing, IR+-valued process. The inverse local time process, also known as 
the ascending ladder time process, is defined as L~^ = {Lf^,t> 0} where 



The ladder height process H = {Ht,t> 0} is defined as 



if t ^ Lqo , 
otherwise. 


The inverse local time L ^ corresponds to the real times at which new maxima are 
reached, and the ascending ladder height process H corresponds to the set of new maxima. 
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The bivariate process {L ^,H) = {(Lj ^,Ht),t > 0}, called the ascending ladder pro¬ 
cess of X, is a two-dimensional (possibly killed) subordinator with joint Laplace transform 

a,/3 > 0. 


The joint Laplace exponent is given by 

K{a,f3) = k{0,0) + adL + Pdn + f {1 — e~°‘^~^^)A{dx,dy), a,/3>0, (2.9) 

-J (0,cso)^ 

where {dL,dH) S 1R+ and A is a bivariate intensity measure on (0,oo)^ satisfying 

/ {1 A\/x^ + 2 /^)A(dx, dy) < oo. 

dro,oo)2 


When L~^ and H are independent, A takes the form A(da;,d?/) = Ai(da;)5o(d?/) -I- 
A_fr(dy)(5o(da;) for x,y>d. In particular, if X is a spectrally negative Levy process, 
one can choose Lt = Mt, which implies that , Ht = Xj,+ =t on {t < Loo}, and 

further K{a, (3) = <I)(q!) -|- /3. 

By letting /3 = 0 in (2.9), we obtain the Laplace exponent of the ascending ladder time 
process, 


-ilog]E[e l{t<i^}] = K(Q;,0) = K(0,0)-|-Q;dL-|- / (1 - e ““)i/L(da:), a > 0, 

^ J(0,oo) 

where VLidx) = A(da;, (0, oo)) is the jump measure of L~^. It follows from integration by 
parts that 


K(a, 0) 


k(0,0) 


= a 




e-°^^VL[x)dx 


a > 0, 


( 2 . 10 ) 


where vl{x) := vl(x, oo). 

The renewal function h associated with the ladder height process H is defined as 

h{x) := [ F{Ht<x}dt, x>0. (2-11) 

J (0,cso) 

When X is a spectrally negative Levy process, it is easily seen that h{x) = dt 

for a; > 0. We recall the follow results (see Theorem 5 in Chapter III and Theorem 19 
in Chapter VI of Bertoin [4]) on the connection between the renewal function and the 
creeping property. Here we say X creeps across x if it enters (a;,oo) continuously. 

Lemma 2.2. The following assertions are eguivalent. 

(i) P{X creeps across x} > 0 for some x > 0. 

(ii) The drift coefficient dn > 0. 
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(iii) The renewal function h is absolute continuous and h' is bounded. 

Moreover, when these assertions hold, there is a version h' that is continuous and 
positive on (0,oo). Finally, lima^^o^ >0 and P{X creeps across a;} =dHh'{x) 
for all X > 0. 

3. Magnitude of drawdowns revisited 

In this section, we revisit some known results of the magnitude of drawdowns of spec¬ 
trally negative Levy processes via fluctuation identities and an approximation approach 
introduced by Lehoczky [24]. Such approach is in the spirit of the general Ito’s excursion 
theory. 

Lemma 3.1. For q>0 and a: > 0, we have 


E[e 


-QTf 


L{m,„>x}] =exp| 




(3.1) 


Proof. For fixed a: > 0 and n G N, let {sn,i,i = 0,... ,n} be a sequence of increasing 
partitions of the interval [0,a;] with 0 = Sn,o < Sn,i < • • • < = x and such that A„ = 

maxi<i<„(s„^i — Sn,i-i) decreases to 0 as n —>■ oo. Using the strong Markov property of X, 
we propose to approximate the event (M,-^ > x) by i 

and thus use 

n 

L'ra:=]^E[e [Aq = s„y_i], 

i=i 

as an approximation of E[e“®L^l{MT„>a:}]- By ( 2 . 4 ), we have 


K=n 


W^^\a) 




= exp-i 




■ + Sn,i - Sn,i-l) - (g) 

VF(«)(a-|-S„y -Sn,i-l) 


Since G 17^(0,00) and is increasing on (0,oo), we have 

f lU(^)(a + - s„.._i) - 1U(«) ia)\\f 1U(«)(a + A„) - (a)\^ , ,2 

IU(^)(a + s„y-s„.._i) J -{ W(«)(a) 

for all 1 < i < n and some constant K > 0. By the fact that — ln(l — e) = £ -|- o(e) for 
small e > 0, it follows that 


E[e l{M,„>x}] 


W^ ‘>\a + Sn,^-Sn,^-l)-W^‘>')ia) 
W^‘1') {a + Sn,i - Sn,i-l) 
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= lim exp s — 7 

n.—^CyO 1 f ^ 


+ Sn.* - Sn.*-l) - 




VF(«)'(a) 

= exp< —TTTTTT^^^ 


W(^\a) 

which completes the proof. 


□ 


By letting g = 0 in (3.1), it is easy to see that follows an exponential distribution 
with mean W{a)/W'{a). Then it follows from (3.1) that, for g > 0 and a; > 0, 


E[e =x\= exp| — 


/lT(«)'(a) W'{a) 


VVF(«Ha) W{a) 

Next, we consider the following lemma which relates to downward exiting. 

Lemma 3.2. For q,s>0, we have 

-,T--s(a-x^.) W{a) zi^\a)wi^^\a)-pW^^\af 

]Ea[e ^0 \Tq <T+] = . 


(3.2) 


IT'(a) 


W^P\a) 


(3.3) 


where p = q — ')p{s), We^'^ and Z'f^ are p-scale functions under P®. 


(p) 


Proof. We first consider that s < ‘h(g), or equivalently, q > ip{s). For 0 < x <y, since 
Tq Ar+ is a.s. finite, by change of measure (2.3) and (2.5), 




-9^0 -s{x-X ) 


° 1 {T”<T+}] - 1 {T-<T+}] 


= ZiP\x) - ZiP\y) 


Ws^P\x) 


(3.4) 


It follows from (2.4) and (3.4) that 

^^^^.,r--sia-x^.)\T_ <T+^ 


= limEaie 

e4.0 


-qTg -s(o-X^_)|^_ 


v'|ro-<T++j 


■.\im(ziP\a)-ZiP\a + e) 

EloV W^^\a + e)J Wia + e)-W{a) 


W{a) ziP\a)W^P^'{a)-pW^P\a)^ 
W'{a) w^^\a) 


(3.5) 
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The other side of the approximation lim£^olEa-e[e ° ^ ^|Tq < T+] also results 

in (3.5). The proof is then completed through an analytical extension of (3.3) to s > 0. □ 

To obtain the main result of this section, we notice that a sample path of X until Ta 
can be splitted into two parts: the rising part and the subsequent crashing part. Because 
of the regularity of 0 for (0,oo), we know that the last passage time = x) = 

{T+\Mr. = x), P-a.s. (see also discussions on page 158 of Kyprianou [20]). Our analysis 
essentially follows this idea: relations (3.2) and (3.3) correspond to the rising and the 
crashing part, respectively. The following quadruple LT is obtained by pasting these two 
parts at the turning point . 


Theorem 3.1. For q,r,s,5>0, we have 

WG+^\a) 

51T0+’')(a) + lT0+0'(a) 

where p = q — ip{s). 


Ej'g-'J'raj 


W^P\a) 


M! 

(3.6) 


Proof. By conditioning on the event (M^^ = x) for x > 0, we have Ta = Gr^ + T~_^ o 
and T~_^ o < T+ o , P-a.s. where 9 is the Markov shift operator defined as 
Xto 9s= Xt+s- Therefore, by (3.2) and (3.3), we obtain 

= \Mr, = x] 

= Vx'-a o < T+ o = x(3.7) 

= V.'-a < T+] 

f /lT0+0/(a) lT'(a)\ ] W(a) Z^P\a)Ws^P^'(a) - 

'vW^0+0(a) ■ lT(a)r/B^ IT?) (a) ’ 

Multiplying (3.7) by the density of Mr^ and then integrating with respect to x, we 
obtain (3.6). □ 


Relation (3.6) generalizes Theorem 1 of Avram et al. [1] by incorporating the joint 
LT of Gr„ and . Moreover, by a similar approximation argument, one can solve for 
the joint distribution of equation (3.6) but with the law of Yt^ , which then recovers the 
sextuple law in Theorem 1 of Mijatovic and Pistorius [29] (the running minimum at 
can also be easily incorporated). 






On magnitude, asymptotics and duration of drawdowns for Levy models 

4. Asymptotics of magnitude of drawdowns 


11 


In this section, we investigate the asymptotics of the LT (3.6) of the magnitude of draw¬ 
downs as a I 0 for spectrally negative Levy processes. Furthermore, we show that such 
asymptotics are robust with respect to the perturbation by positive compound Poisson 
jumps. 


4.1. Spectrally negative Levy processes 

Henceforth, we make the following assumption on the behavior of the scale function at 

Q+. 

Assumption f.l. 

WmxW ix) = 0. 

a:iO 


In fact, since xW’[x) > 0 for all a; > 0, as long as W is well-behaved at 0-1- in the sense 
that 

lima;IF'(a;) = c for some c £ [0, oo], 

a;4.0 

one deduces from the integrability of W at O-f that c = 0. 

Remark f.l. From Lemma 2.1, it is clear that Assumption 4.1 holds if the Gaussian 
component ct > 0 or n(— oo,0) < oo. Moreover, the spectrally negative a-stable process 
with index a £ (1,2), whose Laplace exponent '4>{s) = s“ and scale function 

2 .“— 1 

W{x) = l{a;>0} j 


also satisfies Assumption 4.1. 


Since scale functions are only known in a few cases, we examine sufficient conditions 
on the Laplace exponent to identify cases when Assumption 4.1 holds. 


Remark f.2. For a general spectrally negative Levy process with Laplace exponent 
if, by Lemma 2.1, one can choose an arbitrary sq > 4>(0) and define a function g{x) := 
l{ 2 ,>Q}e“®”^a;IF'(a;), which is non-negative and continuous on ]R\ {0}. By Lemma 3.3 of 
Kuznetsov [18] and (2.6), we further know that g{x) £ L^(R). By integration by parts 
and analytical continuation, one obtains that 



dx = (/7(so ~ is), 


s £ K, 
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where ^{s) := Re(s) > 0. By the Fourier inversion and the dominated 

convergence theorem, we know that a sufficient condition for Assumption 4.1 to hold is 
that (/j(so ~ *•) G A^(K.) as it implies that g{-) is continuous over R. 

Lemma 4.1. Under Assumption 4-1, we have \\mx\a ix) = 0 for every q>0. 


Proof. Since the scale function W is supported on (0,oo), for any k>l, we have 


dx 


W*ik+i)(^x)= [ W'ix-y)W*'^{y)dy + W{0+)W*^ix) 

W'{x-y)dy + W{0+)^ 


< 


x’^-^ 

Jk^. 


W'^+^x), 


(4.1) 


where the inequality above is due to equation (8.23) of Kyprianou [20] and the mono¬ 
tonicity of W. By (4.1) and taking derivatives term by term to the well-known identity 
W^'^'^(x) = where is the fcth convolution of W with itself, we 

obtain 


OO J 

xW<'^'>'(x) = xW'{x) + x'^q'^-^W*^'^+^\x) 

k^l ^ 

= xW'{x) + qxW^{x)e‘>^^^^'>. 


This ends the proof as the right-hand side of the last equation approaches 0 as x 0 by 
Assumption 4.1. □ 


Lemma 4.1 is paramount to derive the following asymptotic results. 

Theorem 4.1. Consider a spectrally negative Levy process X satisfying Assumption 4-U 
For any q,s>0, we have 

if X has unbounded variation, 
if X has bounded variation. 


VF(«)'(e) 
lim — , ,, , 
64.0 VF(9)(e) 


(l-]E[e 


-q'Te-sY.r, 


]) = 


S, 


s + 


q-if{s) 


Proof. Using (3.6), one deduces that 


H^(<?)'(e) 

lU('?)(e) 


(l_E[e-«^-^-'*^--]) 
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= s - o - ' 


+ s{q-f;{s)) [ e-^^W<''^\x)dx + iq-if{s))e-^^W<''^\e). 

J{0,e) 

From the monotonicity of we have 

It follows from (4.2) and Lemma 4.1 that 

Um -^^7^(1 - = s+{q- fj{s))W^‘>^ (0+), 


(4.2) 


which ends the proof by Lemma 2.1. 


□ 


4.2. A class of Levy models with two-sided jumps 

Next, we consider a class of Levy process with two-sided jumps of the form 

Xt=Xt + S+, (4.3) 

where X a spectrally negative Levy process satisfying Assumption 4.1, and 5'“'' is a 
compound Poisson process with arrival rate A+ = 11(0, oo) G (0,oo) and i.i.d. positive 
jump size with distribution function F~^. The two processes X and 5'“'' are assumed to 
be independent. Since we assume that |A| is not a subordinator and is regular for (0, oo), 
it is clear that the same holds for X. 

The characteristic exponent of X is given by 

pOC 

d>{s) = -^{s) + X+ (l-e*'’“)F+(da;), s G R, (4.4) 

Jo 

where ^(•) is the characteristic exponent of X. Henceforth, we add the symbol ~ to all 
quantities when they relate to the spectrally negative Levy component X only. 

By conditioning on the first positive jump arrival time and the jump size, we have the 
following representation of the joint Laplace transform of (Te,!!^^)- 


Lemma 4.2. For q,s>0 and e > 0, we have 


E[< 


,-(JTe-Sy.rEl _ 




{Ts>«;^.Ji+<n+ }J 


1 - (A+/(<z + A+))(l - E[e-(«+^+)-=]) + E[e-?«^+ }] ’ 

(4.5) 

where and are the time and size of the first upward jump of X, respectively. 
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Proof. Recall that is exponentially distributed with mean 1/A+. By the strong 
Markov property of X and the fact that (t^ < < ^t) a-S-, 




{r.>?+,J+>V+_}J 


+ E[e } 


>=] +E[e-^«^ 




{re>«^J+<y+ }J- 
~ 


Solving for E[e one obtains 

E[e-('3+^+)'fe-sfv,] 


E[, 


,-qTs-sY^^] _ 


{re>J+,J+<F+_}J 


1 E[e «€i l{f^>{+,j+>F , _}] 


For the denominator on the right-hand side of (4.6), we notice that 




{fe>?r.jr>V }J 
^ - 


= - E|e-«C 1|,,^ *,] - E|e-«C 1,,., ,] 


-(1-E[e 


q + X+ 

The proof of (4.5) is completed by substituting (4.7) to (4.6). 


(4.6) 


(4.7) 


□ 


We present an analogue of Theorem 4.1 for the Levy process (4.3) with two-sided 
jumps. Note that by (4.4), the drift of the characteristic exponent d of X and X are the 
same when X has bounded variation. 

Theorem 4.2. Consider the Levy model (4-3)- For g, s > 0, we have 


™ W(9+^+)(e) 


(1-E[, 




]) = 


, if X has unbounded variation, 

Q — 

s H---, if X has bounded variation. 


Proof. Since X and 5'+ are independent and (r^ < ^j^) = (fg < ^j^) a.s., we have 

P{r, > ftJt < V-} = 
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<P{f,>^+,J+<e} 

= (l-E[e-^^^'])P{J+<e} 


(4.8) 


It follows from Theorem 4.1 that 


for small e > 0. By (4.5), (4.8) and (4.9), one obtains that 

V7(9+A+)/(^) 


(4.9) 


1(1-E[, 


,-qTe-sY^, 


]) 


]T(9+^+)(e) 

^) 1^^(1 - ^(1 - E [ e -(^+^+)-^])) 


(4.10) 


^(1) 


l-^(l-E[e-(^+^+H)+o(l) 


We hrst consider X has unbounded variation. From Lemma 2.1, we deduce that 
—>■ oo as e), 0. By Theorem 4.1, this further implies that 


1 — E[e 


-(q+A+)fel _ 


] = o(l). 


(4.11) 


One concludes from (4.8) and (4.11) that the denominator on the right-hand side of 
(4.10) approaches 1 as £4,0. Moreover, by Theorem 4.1, 


e™ VF(9+^+)(£) 


(l-E[e-«^=-"^--])=hm 


^(9+A+)/(g) 


eio ITO+^+)(e) 




(1-E[, 


]) = s. 


When X has bounded variation but the Levy measure n(— oo, 0) = oo, note that (4.11) 
still holds by Lemma 2.1. Hence, it follows from (4.8) that the denominator on the right- 
hand side of (4.10) also approaches 1 as e 4,0. Furthermore, by Theorem 4.1, we obtain 


elo ■*')(£) 


lim 


(1-E[, 


,-1 -re_ 


-=]) = s + 


g-b A+ - ip(s) 


X+ g + A+ q-ipis) 

- =S-|-;-. 


q + X+ 


Finally, when X has bounded variation and n(— oo,0) < oo, by Lemma 2.1 and Theo¬ 
rem 4.1, 

















16 


D. Landriault, B. Li and H. Zhang 


Then, by (4.10), (4.8) and Theorem 4.1, it is straightforward to verify that 


™ lF(9+^+)(e) 
which completes the proof. 


(1-E[, 




"=]) =S + 


g-^(s) 


□ 


5. Duration of drawdowns 


In this section, we examine the duration of drawdowns via the LT of the stopping time 
rjt defined in (1.1) for the Levy model with two-sided jumps (4.3). 

To do so, we use the perturbation approach which has been developed by many re¬ 
searchers in similar contexts (e.g., Dassios and Wu [12], Landriault et al. [23], Li and 
Zhou [25], Loeffen et al. [26], and Zhang [37]). To present the main idea, let e > 0 and 
define the following sequence of stopping times: 


=Tr 


’^l = rl+T+ 


o 9 t 


, T* = •d), -f Te O 6», 




S’- — X® 
^0 — 


. 7^+ 


O 0 . 


for f £ N where we recall 0 stands for the Markov shift operator. An approximation of rji, 
is given by 

r]l = inf{t G (r*, iIq] : t — Tg > 6 for some i G N}, 

for which only excursions of Y with height over e are considered. By construction, it is 
clear that Tyf is monotonically decreasing as e 0 , and rjb = lirng^o P-a.s. 

For hxed q> 0, we consider an independent exponential rv with mean 1/q. By the 
strong Markov property of X, 


P{eg >Vb} = lP{eg > > tI -f 6 } -f P{e, > 77 ^,+ b} 

= P{eq A > tI +b} + < e, A + 6 )}P{eq > gl}, 


which yields 


P{eg >b 6 } 


P{eq A'dl> Tg -I- b} 

1 - P{i9j) < e, A (t 1-f 6)} ■ 


(5.1) 


By conditioning on Y^i and then using the strong Markov property of X at time t, 3, we 
find 


PK<e,A(ri+6)} 


[ E[e-«"M{y^^gd,}]P{T+<e,A6} 

J [£, 00 ) 

E[e-«"=] - [ ed,}]P{r+ > e, A 5} 

J [£, 00 ) 


(5.2) 
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and 


P{e, Adl>T^+b}= [ E[e-«"n{K.,6dy}]P{e, A T+ > b} 

J [£, 00 ) 

= 6-^" [ E[e-^-n{^^^ed,}]P{T+>6}. 

J [£, 00 ) 

Substituting (5.2) and (5.3) into (5.1), we obtain 


(5.3) 


E[e-«’'^']=PK>%"} 


1 - E[e-9-=] + E[e-9-. l^y^^ gdy}]P{r+ > e, A 6} ' 


(5.4) 


From the representation (5.4), it seems relevant to define, for a; > 0 and p > 0, a bounded 
auxiliary function 


/4(t) := [ E[e-^"n{y^^edy}]P{T+ > A t} 

J [£,00) 

= f E[e"'?^'l{K,^gdy}]P{A7epAi<2/}, 
J [e.oo') 


(5.5) 


where the dependence of (5.5) on q is silently assumed. Hence, we rewrite (5.4) as 


= ^ —. (5.6) 

l_E[e-9-.] + /4(6) 

To obtain a well-defined asymptotics for fe^\b) as e4.0, the key is to investigate the 
convergence of the measure P[e“'^’^''l{Y^^gdi/}] as e ), 0, which is closely related to the 
asymptotic results of Section 4. As we will see below, the convergence of the measure 
differs according to whether the Levy process has bounded or unbounded variation. 


5.1. Bounded variation case 

We hrst show that the distribution function of the running maximum of X is well- 
behaved. 

Proposition 5.1. Let X be a Levy process of bounded variation with a drift d> 0 in 
its characteristic exponent representation (2.8). Then, for any fixed p>0 and t>0, the 
function PjMepAt ^2/}/?/ bounded for y G (0,oo). Moreover, if we further assume that 
n(— oo,0) = 00 and H has no atoms on (—oo,0), the function ¥{Me^/\t < y}/y is also 
continuous for every y G (0, 00 ). 
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Proof. For any fixed p>0 and t > 0, we denote by 

F^^P\y) :=P{Me,At<2/}/2/. 


We first consider the case p = 0. Using the upper bound in equation (4.16) of Chaumont 
and Malecki [7] (which holds for a general Levy process), we know that 


Ft°\y) < 




/ y 


(5.7) 


where we recall h{-) is the renewal function defined in (2.11). Since X has bounded 
variation and d > 0, we deduce that X creeps upwards by Theorem 7.11 of Kyprianou 
[20]. From Lemma 2.2 we know that h{y)/y converges to a finite limit as y 0. Therefore, 
we conclude from (5.7) that F^^\y) is bounded for y G (0,oo). 

Next, we consider the case p > 0. By Wiener-Hopf factorization, it is well known that 
Mep follows an exponential distribution with mean l/<l>(p) >0. Moreover, since Mt > Mt 
a.s. for any t>0, one obtains that 


J(o,t) y y 

f p^-ps ^{Ms<y} 

J fo,oo) y 

.e-^‘uf)(y). 


'(o.t) 

/ 

' (0,oo) 

1 _ g-i(p)y 


By the boundedness of F^^\-), we deduce that F^^\y) is also bounded for y G (0,oo). 

Finally, suppose that we also have n(— c»,0) = oo and 11 has no atoms on (—oo,0). 
For any fixed t > 0, by Theorem 27.7 of Sato [33], we know that the law of Xt is absolute 
continuous with respect to the Lebesgue measure, so is the law of Xt from the property 
of convolutions. In addition, by Theorem 6.5 of Kyprianou [20], we know X is regular for 
(0,oo) as X has bounded variation and d > 0. Therefore, from Theorem 1 of Chaumont 
[6], we conclude the law of Mt is absolute continuous with respect to the Lebesgue 
measure. As a consequence, P{MepAt ^y}/y is continuous for every y G (0,oo). □ 


Remark 5.1. For the Levy model (4.3) with X has bounded variation and n(— oo, 0) = 
oo, it follows that P{MepAt < y}/y is bounded and continuous for y G (0,oo) due to our 
assumptions that |A| is not a subordinator, X is regular for (0,oo), and B has no atom 
on (— 00 , 0 ). 


We are now ready to present the main result of this subsection. 
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Theorem 5.1. Consider the Levy model (4-3). If X has bounded variation and satisfies 
Assumption 4-1, for any q> 0, we have 




1(0,00) ^{Mb<y}U{-dy) 


9 + /(o.oo) IP{^e,Ab < 2/}n(-d2/) ■ 

Proof. We first consider the case n(— oo,0) = oo. From (5.5) with p > 0, we have 


^r{q+X+),^^^ 


f(p\b) 


^(9+A+)/(£) 


W(9+^+)(£) 7[,.oo) 


/ ]E[e '^^"l{K,^6dy}]IP{MepAf)<2/} 
J k.oo) 


' (0,oo) 


P{Me,Ab<?/} , 

l-e-v ■ W('?+^+)(e) ^ Im.Gdy}] 


(5.8) 


' (0,oo) 


PjM^BpAb < y} 
1 - e-y 


h'eidy), 


where p-sidy) is a finite measure on (0,oo) defined as 


Wi<!+X+)'(^) 

h-eidy) = |^(q+A+)('g.^ ® ^)lE[e 'l{Kr^Gdy}]- 


(5.9) 


By Theorem 4.2, we have 


r y^(q+\+)l l^\ 

_ 1^(s) - ^^(s + 1) 

^ d 


for all s > 0. On the other hand, we notice from (2.2) that 



1 - e~y 

e-^y——Ili-dy) 

d 


1 

d 


' (— 00 , 0 ) 


midy) - i 

“ J{-oo,0) 


l)Il{dy) 


^ ^ ^(s)-^(s + l) 
d 


Hence, by Proposition A.l, one concludes that, as e4.0, ^e(dy) weakly converges to the 
measure d“^(l — e“*')n(—dy), which is a finite measure on (0, oo) because X has bounded 
variation. 
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From Proposition 5.1 and Remark 5.1, we know the function PjMepAt) < 2/}/(l ~ is 
bounded and continuous for y € (0,c»). By the definition of weak convergence, it follows 
from (5.8) that 


p(y(g+A+)// p 

lim— rTTTTT^f^^Hb) = lim / 
eio 1F(9+A+) (£)■'" ^0 7(0, 

-L 


1-e-y 




l-e-y d^ ^ ^ 


I 


(0,oo) 


IP{A^epA6 < 2/}n(-dy). 


Therefore, by (5.6), (5.10) and Theorem 4.2, we have 


E[e-™1 = _ e;0 () _ 

W(.+x+)(e) (l - w(-^+>+)(e) ^ (^) 


e I^o,oo)^{^b<y}Ti{-dy) 

9 + /(O.oo) IPi^e^Ab < ?/}n(-d2/) ■ 


(5.10) 


Finally, we consider the case that n(— oo,0) < oo. By (4.12), for any s > 0, we 


have 


lim f e "^E[e 

J(0,oo) 


'0(s) — sc? + n(— c», 0) 
q + n(— oo, 0) 


f n(-di,) 

V(o,™) 9 + n(-co, 0 ) 


By Proposition A.l, we see that the measure l^y^^^dy}] weakly converges to the 

measure n(— dy)/(g + n(— oo, 0)) as e ), 0. Since PjMepAfc < 2/} is bounded and upper 
semi-continuous in y G (0, oo), it follows from Portemanteau theorem of weak convergence 
that 


limsup/JPi(6) = limsup /" P{MepAb < y}]E[e 

£■4-0 S-i-O J (0,oo) 


< 


1 

q + n(— oo,0) 



P{AfepAb < 2/}n(-dy). 


(5.11) 


On the other hand, since P{Mej,Ab < u} is lower semi-continuous in y £ (0,oo). By Porte¬ 
manteau theorem again, we have 


lim inf f'f\b) > lim inf [ P{Me^Ab < y}E[e gdy}] 

£4,0 £40 7(0,00) 

]P{AfepAb < y}U{-dy) 


g-bn(-oo,0) 7(0,oo) 


(5.12) 
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1 

q + n(-oo,0) 



P{MepAb < 2/}n(-dy), 


where the last equality holds because 11 has no atom on (—oo,0) and PjMepAh < y} = 
P{MepAb < y} for almost all y > 0. By letting e 10 in each term of (5.6) and using (5.11), 
(5.12) and (4.12), this completes the proof of Theorem 5.1. □ 


5.2. Unbounded variation case 


We now consider the unbounded variation case for which the following assumption on 
the density of Xt is made. 

Assumption 5.1. If X has unbounded variation, we assume that the density of Xt, 
namely p^{x), is hounded for all t > 0. 

Remark 5.2. We point out that Assumption 5.1 is identical to assumption (HI) in 
Chaumont and Malecki [7], which is equivalent to the assumption that the characteristic 
function G L^(R), for all t>0. It is also clear that, if A is a spectrally negative 

Levy process with unbounded variation and Y is an arbitrary Levy process independent 
of X, then the sum X + Y satisfies Assumption 5.1 as long as X does. Hence, examples 
of Levy processes satisfying Assumption 5.1 include processes with cr > 0, or cr = 0 and 
with a spectrally negative a-stable jump distribution with a G (1,2). 

The following proposition shows that, for a Levy process with unbounded variation 
satisfying Assumption 5.1, the density of the running maximum at 0+ is well-behaved. 


Proposition 5.2. Let X he a Levy process with unbounded variation that creeps upwards 
and satisfies Assumption 5.1. Then the running maximum Mt has a continuous density 
p¥{') fo'f' every t > 0 and further. 


limpf(x) 

a; 4-0 


-I- k(0,0) 

dn 


> 0 , 


where vl{‘) is the tail of the jump measure of the ascending ladder time process (see 

( 2 . 10 )). 


Proof. From Lemma 2.2, we know that the renewal density h' can be chosen to be 
a continuous function with well-defined limit h'(0) = ^ > 0. Since X has unbounded 
variation, assumption (H2) of Chaumont and Malecki [7] also holds. By Proposition 2 
and Theorem 1 of Chaumont and Malecki [7] , we know that Mt has a continuous density 
pf^{x) for every t > 0, and also, 

(^) = ^(C > 0, 

xiO n yx) iciO 
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where n((^ > t) is the excursion measure of excursions with length over t>0. From (6.11) 
and (6.14) of Kyprianou [20] (see also Section IV.4 of Bertoin [4]), we know that 

n(C >t) = VL{t) + k(0, 0) > 0, 


which completes the proof. 


□ 


Corollary 5.1. Under the eonditions of Proposition 5.2, for any fixed p>0 and t > 0, 
the function PjMepAt < y}/y *5 bounded and continuous for y € [0,oo), where its value 
at y = 0 is defined as the right limit 

F{Me /\t<y} ^ f f -PS- / , -pt- 

lim-= — / pe ^ VL(s)(ls + e ^ j^L(t) + k(0,0) . 

2/10 y dn \J(o,t) ) 


Proof. From Proposition 5.2, it is only left to justify the limit of PjMepAt ^y}/y as 
yfO. By dominated convergence theorem and Proposition 5.2 again, we have 


lim 


n^^epAt < y} 

y 


f 




pe lim 


■ ds + e ^ lim ■ 

yfo y 


pe limp^( 2 /) ds + e ^‘hmpf^(y) 

viO ylO 


pe PVi(s)ds + e ^^^(t) + k(0,0)^ 


which ends the proof. 


□ 


Now we are ready to present the main result of this subsection. 

Theorem 5.2. Consider the Levy model (4-3). If X has unbounded variation and sat¬ 
isfies Assumptions 4-b 5.1, for any g > 0, we have 

Tpfp-OT!.] _ p-qb _ ^Ljd) + s^(0,0) _ 

1(0,b) qa-i^uLft) At + e-^’^VLib) + k(0, 0)' 

Proof. It is clear that the Levy model (4.3) creeps upward as its spectrally negative 
component X does and its upward jumps follow a compound Poisson structure. More¬ 
over, since X satisfies Assumption 5.1, by Remark 5.2, we see that all the conditions of 
Proposition 5.2 are satisfied. 

For the finite measure /ie(d?/) defined in (5.9), it is straightforward to verify from 
Theorem 4.2 that, for any s > 0, 
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= 1= [ e-^ySo{dy). 

JR+ 

It follows from Proposition A.l that Heidy) weakly converges to the Dirac measure i5o(dy) 
as £4^ 0- Moreover, by Corollary 5.1, we know that the function F{Me At < y}/(l — e“^) 
is also bounded and continuous for y € [0,oo), where its value at y = 0 is defined by the 
limit as yJ-O. From (5.8) and Corollary 5.1, we have 


eio lV(9+^+)(e) eloL 1 - e-v 


= lim 


^{MepAb < y} 

1 - e-y 


(5.13) 


dn \J{o,b) 

It follows from (5.6), (5.13) and Theorem 4.2 that 


^( [ Pe y*i^L{t)dt + e P’^VLib) + k{0,0)\ 

'-H \J(0,b) ) 




lim, 


w 


f(0) 




fr(b) 


= e-y» 


limeio (1 ~ + fe'^hb)) 

_ VLib) + k(0,0)_ 

f(o,b) ye-‘?*pL(t) dt + e-y^pL{b) + k(0, 0) ’ 


which ends the proof. 


□ 


In general, the function and k( 0, 0) are only implicitly known via (2.10) and Wiener- 
Hopf factorization. When X has no positive jumps, we can express explicitly in 

terms of pf. 


Corollary 5.2. Let X be a spectrally negative Levy process with unbounded variation 
and satisfies Assumptions 4-1 rtnd 5.1. For any q > 0, we have 


= e"«'’ 


/(f,,oo) ipf(0)ds 


/( 


( 0 .&) 


qe 


/(t,oo) iPf (0) ds dt + e-9^ ipf (0) ds ■ 


Proof. By Kendall’s identity, for any fixed t,y >0, we have 


-F{Mt <y} = - [ P{T+ e ds} = [ -yf (y) ds. 

y y J{t,oo) J{t,Qo) ^ 


It follows from Fourier inversion that, for any y € M and s > 0, 
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From the proof of Proposition 5 of Chaumont and Malecki [7], we know that for any 
fixed t > 0, 




1 


|e dwds < oo. 


(t,G 


By the dominated convergence theorem and Corollary 5.1, we have 

1 


+ k(0,0) 

- - -= lim- 

dn yio y 

which completes the proof by using Theorem 5.2. 


(tjOo) 


Ps (0)ds, 


□ 


6. Examples 

Example 6.1. Consider a Brownian motion, i.e. Xt = pt + aWt with cr > 0. For any 
fixed t > 0, we have 


Pt (x) = 


1 


2 aH 

By Remarks 4.1, 5.2 and Corollary 5.2, we have 


[x — ytY 


E[e 




1(0,b) ge-i*g(t)dt + e-i’>g(b)’ 


where g(t) := jp^(0)ds = ^ ^ */(2o- ) _ and N(-) is the cumula¬ 

tive distribution function of a standard normal rv. 


Example 6.2. Consider a spectrally negative a-stable process with Laplace exponent 
ip{s) = s“ with a G (1,2). For fixed t > 0, it is well known (e.g., pages 87-88 of Sato [33]) 
that 

pf (^) = sinf—) {t-^/-xr-\ 

7 T ■“ n\ \ Oi J 


n—l 


where r(-) is the Gamma function. It follows that 


f -pf{0)ds = -Tf-']sin[-']t- 

J(t,oo) s tt \aj \aj 


' (t,oo) 

By Remarks 4.1, 5.2 and Corollary 5.2, we have 


E[e 


-mb] — 


1 


gqrfp^l/a gg-gt^-l/a dt 1 ' 
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Example 6.3. Consider a spectrally negative Gamma process with Laplace exponent 
'0(s) = sd+ f — l)/3|a;|“^e““ dx = sd —/31og(l + s/a), sSHI''', 

J ( —cso,0) 

where a,/3 > 0 are constants. From Remark 4.2, we define 

, ,_s^'(s) - if{s) _ /?log(l + s/a) - /?s/(s + a) 

' 'il){sY (sd —/31og(l + s/a))^ 

One can easily verify that, for any fixed sq > 4)(0), we have p{so + i) € L^(]R) which 
implies Assumption 4.1 holds. Using Kendall’s identity and 

Pti.x) = xGR,t>0, 

we have 

^{Ty>t} = yf -■p^^l{y<sd}(srf-d)^""^e"“(®‘'"^)ds, y > 0 and t > 0. 
d(t.oo) s r(/3s) 

Hence, by Fubini’s theorem followed by some calculations, we can show that 


(0,oo) 


^{Mt < y}n(-d?/) 


[ F{T+ > t}U{-dy) 

J (0,oo) 

f f ^ p~^y 

/ y -TTr;TT^{v<sd}{sd-yf^-^e-°‘^^‘‘-y''ds(3 - dy 

J( 0 ,oo) J(t,oo) sT{f3s) y 


f (0,oo) */ (i,oo) 

= f 

/(t.oo) r(^s) 
Using Theorem 5.1, one concludes 


3s-2 


E[i 




^fis-2^-a,sd 

T( 0 S) 


ds 


q + (da)/5« /(p dt ds + e-^\daY^ ds' 

Example G.f. Consider Kou’s jump-diffusion model given by 


n: 


n: 


Xt — pt + aWt + ^ ^ Jj 


i=i j=i 
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where n gM., a > 0, are two independent Poisson processes with arrival rates > 0 
and are a sequence of i.i.d. exponentially distributed random variables with mean 
1 / 77 ^ > 0. Its Laplace exponent is given by 

Hs) = y + MS + A- - 1) + A+ - 1), s € (-77-,^+). 


According to Corollary 1 of Asmussen et al. [1] and Section 6.5.4 of Kyprianou [20], it is 
known that the Laplace exponent of the ascending ladder height is given by 


! Q\ (/5 + Pl,a)(/3 + P2,a) 

=- ,r, , - 

(^ + 77 + 


a, /3 > 0, 


where p\^a and p 2 ,a (with < 77 + < p 2 ,ct) are the two distinct nonnegative solutions 
of '0(s) = a. By Remarks 4.1, 5.2 and Theorem 5.2, one obtains 


E[e 


-9»?6l _ g-9b 




/(o.b)9e '?*7/L(t)dt + e '?Vi,(6) + 


P1,QP2,0 • 






Appendix 

The following result is from Theorem 5.22 of Kallenberg [17]. 

Theorem A.l (Extended continuity theorem). Let /ii,M 2 )--- probability mea¬ 
sures on with characteristic functions fin{t) —>■ ip{t) pointwisely for every t G where 
the limit p is continuous at 0. Then /in converges weakly to p, for some probability mea¬ 
sure p in with p = p. A eorresponding statement holds for the Laplace transforms of 
measures on . 

Proposition A.l. Let {/injneN be finite measures on [0,oo) with Laplace transforms 

Pn{s)= f e“®*'/i„(dy), 

*/ IR-[- 

for n € N and s > 0. Suppose that lim„_>oo An(s) = 7’(s) for all s > 0, where (p{-) is a 
positive and continuous function on [0,oo). Then pn weakly converges to p as n^oo, 
for some finite measure p on [ 0 ,oo), and p = ip. 

Proof. Since hm„_>.oo Pn{0) = <m( 0) > 0, we can consider a sequence of probability mea¬ 
sures Vn{dy) := pn{dy)pn{0)~^■ By our assumptions, it is easy to see 

lim 7>„(s) = (m(s)m(0)“\ 

n—¥oo 

which is a continuous function at 0. By Theorem A.l, one concludes that {i^nlneN weakly 
converges to some probability measure v on [0, 00 ) with v(-) = p{-)p{0)~^. Therefore, by 
letting /i(-) := i;(-)(M( 0 ), we can see that pn weakly converges to /i as 7 i —>■ 00 and p = p.\3 
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